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Abstract 

We derived the three-dimensional distribution of H I gas in the Milky Way Galaxy 
using the latest H I survey data cubes and rotation curves. The distance of the H I 
gas was determined by the kinematic distance using a rotation curve. We solved the 
near-far problem in the inner Galaxy by a fitting method which involves introducing 
a model of vertical H I distribution. In our resultant maps we could trace three 
prominent arms: the Sagittarius-Carina arm, the Perseus arm, and the Outer arm. 
These three arms were found to be logarithmic spiral arms. The pitch angles of the 
Sagittarius-Carina, Perseus, and Outer arms were estimated to be about 11°, 18°, and 
7°, respectively. The Sun is located in a region rich in H I gas between the Sagittarius- 
Carina arm and the Perseus arm. The H I disk shows large and asymmetric warping 
in the outer disk: the H I disk goes up to about 1.5 kpc above the Galactic plane in 
the northern hemisphere, and down to about 1 kpc in the southern hemisphere, which 
means asymmetric warping. The inner H I disk is also found to be tilting. The radius 
of the H I disk is about 17 kpc and the Hi mass within this radius is estimated to be 
2.5 x 10 9 M , which corresponds to 1.5% of the dynamical mass predicted from the 
rotation curve. We also found that the H I outskirt is largely swelling in the fourth 
quadrant, and hence the Galaxy is significantly lopsided. The scale-height of the H I 
layer increases with the radius, and is correlated with the H I volume density at the 
centroid of the H I layer. 

Key words: Galaxy: disk — Galaxy: kinematics and dynamics — Galaxy: struc- 
ture — ISM: kinematics and dynamics — radio lines: ISM 



1. Introduction 

While the Milky Way Galaxy is the most explored galaxy, its global three-dimensional 
(3D) structure, such as the spiral arms, has not yet been clarified, because we are embedded 
in it. One of the best tracers of the global structure of the Galaxy is neutral hydrogen (H i) 
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gas, because H I gas is one of the major components of the interstellar matter (ISM) in spiral 
galaxies, and the ISM is almost optically thin against the 21-cm line of the H I gas. 

The earliest study of the global structure of the Milky Way using H I data was done 
by Oort, Kerr, and Westerhout (1958), who presented a face-on H I density map of the whole 
Galaxy. They reported several Hi arms in the Milky Way: the so-called Sagittarius arm, the 
Orion arm, and the Perseus arm. Later, as more H I surveys were achieved, Kerr (1969), 
Weaver (1970), and Verschuur (1973) traced larger-radius H I arm, which is the so-called Outer 
arm, on a longitude-velocity diagram. Kulkarni, Blitz, and Heiles (1982) analyzed the H I data 
to make an H I density map for the outer Galaxy, and traced two arms having a pitch angle 
of ~ 22°-25°. Burton and te Lintel Hekkert (1986) showed three-dimensional H I maps of the 
outer Galaxy for all Galactic longitudes, and reported a warped H I disk. Diplas, Savage (1991) 
also presented an H I density distribution of the outer Galaxy, and discussed the warped disk 
and the flaring H I disk. 

More recently, Hartmann and Burton (1997) archived a deeper and more complete survey 
with a much higher sensitivity. In the present work, we tried to convert this latest data cube 
into an H I distribution. Combining with the southern survey data cubes available, we aimed to 
derive a 3D map of the whole Galaxy. This paper is the first of a series presenting our analysis 
of the 3D structure of the ISM. In the second paper, we will present the result for molecular 
gas. 

In this paper we describe the data and our method for deriving the H I distribution of 
the whole Galaxy in sections 2 and 3, and report on the resultant H I maps and some properties 
of H I disk of the Galaxy in section 4. 

2. Data 

2.1. Hi Survey Data 

Three available H I survey data are used in this paper: the Leiden/Dwingeloo survey 
(Hartmann, Burton 1997), Parkes survey (Kerr et al. 1986), and NRAO survey (Burton, Liszt 
1983), which all are archived in the FITS format. We took data within a latitude range of 
—10° < b < +10° from these data cubes to study the H I distribution. 

The Leiden/Dwingeloo survey (Hartmann, Burton 1997) is used for the northern hemi- 
sphere, which was obtained with the 25 m Dwingeloo radio telescope in the Netherlands with a 
half-power beam width (HPBW) of 0.°6. The data consists of 21 cm line spectra over the entire 
sky above declinations of —30° with a grid spacing of 0.°5. The radial velocity (V v ) coverage is 
from —450 to +400 km s _1 with a velocity resolution of ~ 1.03 km s _1 . The r.m.s. noise level is 
0.07 K. These spectra are averaged so that the bin width of latitude becomes 1° in this paper. 

The Parkes survey (Kerr et al. 1986) is used for the southern hemisphere, which was 
obtained with the CISRO 18 m telescope in Australia with an HPBW of 0.°8. It covers Galactic 
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Table 1. Modified coefficients of the inner rotation curve. 





R < 0.09 Rq 


0.09 R <R< 0.45 Rq 


0.45 Ro<R<R 




0. 


325.0912 


-2342.65649 


Ax 


3261.30 


-264.0309 


2663.95435 


A-2 


-17847.8 


261.7602 


-1156.07764 


A, 


52752.5 


-132.822906 


269.354645 


A 4 


-87027.2 


31.9537735 


-36.2039909 


A 5 


74344. 


-2.857958 


2.8025787 


A 6 


-25510. 


0. 


-0.1158279 


A? 


0. 


0. 


0.001977 



longitudes from I = 240° to 350° with a grid spacing of A/ = 0.°5, and Galactic latitude b = — 10° 
to +10°, with a grid spacing Ab = 0?25. The velocity range is 300 km s" 1 with a velocity 
resolution of ~ 2.1 km s _1 . The average r.m.s. is 0.27 K. These spectra are also averaged so 
that the bin width of latitude becomes 1°. 

The Galactic center H I survey of National Radio Astronomy Observatory is used for the 
southern hemisphere between I = 350° and / = 360° (Burton, Liszt 1983). This data cube was 
obtained with the 140- feet telescope of the NRAO, whose HPBW is 21'. The survey covers the 
region 348° < I < 10°, —10° < b < 10°, and a grid spacing of 0?5 in both / and b. The velocity 
range is \v\ < 310 km s _1 and the velocity resolution 5.5 km s _1 . The r.m.s. level is 0.03 K. 
These spectra are also averaged so that the bin width of latitude becomes 1°. 

2.2. Rotation Curve 

Since the 21-cm spectra have only information about the radial velocity of the H I gas, 
we need to transform the radial velocity, V T , of the ISM in the Galaxy into the distance, r, in 
the line of sight. The relation between the radial velocity, V T , and the distance, r, in the line of 
sight is obtained using the rotation curve of the Galaxy, assuming that it rotates circularly. 

For the inner rotation curve, we use Clemens's rotation curve (Clemens 1985), which was 
derived from the Massachusetts-Stoney Brook Galactic equator CO survey. Clemens's rotation 
curve was fitted by a polynomial of the form 

V(R) = XUAR*. (1) 

Though Clemens listed the fitted coefficients, A iy for the cases of (Rq, Vq) = 
(8.5 kpc, 220 km s _1 ) and (10 kpc, 250 km s _1 ), we modified them for the case of (Rq,Vq) = 
(8 kpc, 217 km s -1 ) in order to connect smoothly to the outer rotation curve in this paper. 
Modified coefficients are listed in table 1. 

For the outer rotation curve, we use a model rotation curve of Dehnen and Binney (1998), 
which was determined by fitting to observational constraints ever obtained. However, because 
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of the small number of data points for the outer rotation curve, the mass model has not been 
settled yet, and four types of models are presented by Dehnen and Binney (1998). Of the 
four models, Model 2 is adopted in this work, since it has the smallest deviation from the 
observation among the four. In using this rotation curve, we adopted Rq = 8 kpc, Vq = 217 km 
s _1 , which are used for the Galactic constant in Dehnen and Binney (1998). The shape of the 
whole rotation curve used in this paper is shown in figure 1. 

3. METHOD 

In this section we describe a method for transforming the 21-cm spectra into an H I 
density distribution. Due to a geometrical reason, the H I distribution should be calculated 
separately in several cases: the outer Galaxy, the inner Galaxy, and near the tangential points. 
The radial velocity and the distance in a line of sight can be uniquely determined in the outer 
Galaxy, but the distance cannot be determined uniquely from the radial velocity in the inner 
Galaxy because of the near-far ambiguity. Since the H I gas has a velocity dispersion, the H I 
density at the tangential point cannot be calculated without considering its effects. 

3.1. Basic Idea and Derivation of the Hi Distribution in the Outer Galaxy 

The distance of the H I gas having a radial velocity V T can be calculated as follows. 
Assuming that the motion of the Galactic ISM is purely circular and that the gas away from 
the Galactic plane rotates at the same speed as the gas in the underlying disk, the radial 
velocity V T of the ISM at Galactic longitude I, Galactic latitude b, and Galactocentric distance 
R is expressed by 



V T {l,b) 



^V(R)-V 



sin /cos b. (2) 
The heliocentric distance r of the ISM at (R,l,b) is obtained by solving the equation 

R 2 = r 2 + R 2 -2rR cos/. (3) 



In the outer part of the solar orbit (R > R ), we have one solution, r = R cosl + y R 2 — R^sin 2 l, 
because the other solution has a negative value, which is inadequate. 

In the inner Galaxy (R < Rq) we have two adequate solutions, r = R Q cos I ± 



R 2 - i^sin 2 /; that is, we have the two-fold of kinematic distance in the inner Galaxy. However, 
at a tangential point where R = R^sinl, equation (3) has only one solution despite the inner 
Galaxy. Hence, we have no ambiguity of distance at this point, and the kinematic distance is 
uniquely determined for a given radial velocity. For example, we plot the radial velocity against 
the heliocentric distance (r) in the case of / = 30° in figure 2. The radial velocity in the longitude 
range / < 90° becomes positive in the inner Galaxy, and increases toward the tangential point 
where it gets maximum. Beyond the tangential point it decreases monotonically and becomes 
negative in the outer Galaxy. The radial velocity for the whole Galactic equator is plotted in 
figure 3. 
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The H I density can be estimated as follows. If the 21-cm line is assumed to be optically 
thin, the H I column density TVhi [cm -2 ] is given by integrating brightness temperature, Tb [K], 
over radial velocity, V r [km s -1 ], 

iV H i[cm- 2 ] = 1.82 x 10 18 / T h dV T . (4) 

In this work, the velocity width, AV T , is taken to be 2.06 km s~ x in the case of the 
Leiden/Dwingeloo survey, 2.00 km s _1 in the case of the Parkes survey, and 5.5 km s _1 in 
the case of the NRAO survey based on the velocity resolution of the southern survey. 

The average volume density at an individual point is given by dividing the column density 
by the path length, Ar, which contributes emission in the velocity range V t \ < V r < V?x + AV T : 

n ffl [cm- 3 ] = ^ = 1.82xl0 18 T b ^. (5) 
Ar Ar 

Thus, we obtain the H I densities, rim, i n the outer Galaxy from the 21-cm spectra 
Tb(l,b,V T ). However, the local gas with a high velocity dispersion contributes to the high- 
altitude emissions. Therefore, the obtained H I distribution includes the local gas in the form 
of riobscrvcd = ft-Hi + ft-iocai- I n this paper, we assumed that the local gas can be expressed by 
a linear function of the Galactic latitude, b, and subtracted it from n b se rved- Because of this 
procedure, the diffuse H I gas, which extends to a high Galactic altitude (Diplas, Savage 1991), 
is subtracted, and can not be detected in the maps presented in this paper. 

3.2. Hi Density around the Tangential Points 

At a tangential point r = R cos / the kinematic distance can be uniquely determined 
even for the inner Galaxy. Since the H I gas has an intrinsic velocity dispersion, the H I spectra 
have a higher velocity component than the terminal velocity, V t . Velocity dispersion, a, is 
assumed to be 10 km s _1 ; Malhotra (1995) estimate it to be 9.2 km s _1 for the first quadrant 
(0° < / < 90°), 9.0 km s" 1 for the fourth quadrant (270° < I < 360°). In a given line of sight, the 
H I gas having a radial velocity of the range \ V t \ — a < \ V T \ < \ V t \ contributes to a higher velocity 
component than the terminal velocity due the velocity dispersion. Also, equation (5) at the 
tangential points does not work because AV T /Ar goes to zero there. Therefore, the volume 
density of the H I gas giving the radial velocity |Vt| — c < |V^| < |Vt| around the tangential 
points is calculated by integrating the emission within the velocity range of |Vt| — er < \V r \ < oo, 
as follows: 

, x 1.82 x 10 18 /^ a T h dV r . 

nm r t = -jtt ^ JVt ~" \ , (0° < K 90° , 6 

r 2 {V t - a) - ri(V t - a) 

i.82xio 18 /yg g r b dv r 

nm{r t ) = —— — r — - — -, 270 < / < 360° , 7 

r2{Vt + a)-r 1 (V t + a) 

where ri(Vt±<r) and r2(K±cr) denote the heliocentric distances of the points giving radial 
velocities of V T ± a. Subscripts 1 and 2 indicate the near and far points, respectively. 
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3.3. Derivation of the Hi Distribution in the Inner Galaxy 

In the inner Galaxy there are two points giving the same radial velocity in a given line 
of sight, and thus the observed spectra give the total emission at the near and far points. We 
solved this problem of the distance-degeneracy by considering the H I layer thickness. 

First of all, we determined the H I layer thickness and its variations using the 
Galactocentric distance. For the model of the vertical distribution of H I we use Spitzer's model 
(Spitzer 1942). Spitzer (1942) analytically solved the dynamical ^-distribution of stars and 
ISM. Assuming that the stars and gas in the Galactic disk are isothermal and self-gravitating, 
the number density nni can be expressed as a function of z by solving Poisson's equation in 
hydro-dynamical equilibrium. We then obtain 

"m(0 =n H i sech 2 (£) (8) 
with £ being 

£ = log(l + V2)i-^, (9) 

z l/2 

where z is the height from the Galactic equator, zq is the height where nni becomes the max- 
imum value, riHio) an d Z1/2 is the height where nni becomes half the value of the maximum, 
which is half the full width of the half maximum (FWHM). 

In equation (8) we assume that n-m and z vary with its position, but the scale-height, 
Zi/2, depends on the Galactocentric distance alone. We determined the H I layer thickness 
Z1/2 as a function of the Galactocentric distance by fitting the model function to the observed 
data. The thickness of the inner H I layer was determined by fitting the model function to the 
vertical distribution at the tangential point obtained in the previous section; we substituted 
the values into three parameters, (tibi , Zq, and z 1 / 2 ) in equation (8), and chose the best set 
of parameters to fit to the observation by the least-squares method. The steps to search for 
the parameters were 0.001 of the peak for nm , 1 pc for z and zu2- On the other hand, the 
thickness of the outer disk was determined by fitting the model to all points, since there is no 
distance-ambiguity in the outer Galaxy. The resultant thickness of the H I layer is presented 
in figure 4. 

For the next stage, the total emission was divided into two H I densities at the near and 
far points in the line of sight using the thickness of the H I layer. First, consider an H I gas having 
a radial velocity V < V T < V + AV for a given longitude, I. This gas locates at corresponding 
kinematic distances, r\ — Ar x /2 < r < T\ + Ari/2, and r 2 — Ar 2 /2 < r < r 2 + Ar 2 /2. Assuming 
that the H I gas follows Spitzer's model, the H I column density, Nm(b), is represented as a 
function of latitude with four parameters: riHi 0l , ^hi 02 , z 0i , an d z 02 . That is, 

N m (b) = n m sech 2 (£i) x + n H i sech 2 (£ 2 ) x -, (10) 

1 coso 2 coso 

where 
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a=io g (i+V2) 



ritanfe — z 0l 



(11) 



£ 2 = log(l + \/2) 



Z\/2 

^tanfr — z 02 



(12) 



Zl/2 



The four parameters were determined by fitting the observed data by a model using the 
least-squares method. Finally, the H I density at the near and far points were calculated using 
(^hi 0i ,^0i) and (n H i 02 ,2o 2 ), respectively. 

3.4- Construction of the 3D H I Density Map 

We thus obtain an H I density cube, n m (l,b, r), from the 21-cm spectra T h (l,b, V T ). 
Finally, the nui data cube in the (l,b,r) space was transformed into tzhi in (x,y,z) space where 
x,y,z are the Cartesian coordinates in the Galaxy. The position (x,y,z) = (0,0,0) indicates the 
Galactic center (G.C.), and the Sun is situated at (0 kpc, 8.0 kpc, kpc). The space of (l,b,r) 
was also transformed into cylindrical coordinates (R,8,z), where R and 9 are the distance from 
the rotation axis and the azimuth angle around it, respectively. The angle 9 was taken so that 
9 = 180° points to the Sun and 9 = 90° is parallel to I = 90°. 

4. Results 

The derived H I density distribution is presented in the forms of (1) a face-on H I column 
density map, (2) face-on Hi volume density maps sliced at several altitudes, and (3) edge-on H 
I volume density maps sliced at several azimuthal angles around the Galactic center. 
1. face- on column density map — Figure 5 shows a face-on column density map of H I gas in 
the Galaxy. According to the Spitzer (1942) 's model, the H I column density, S, is calculated 
as 



2. face-on volume density maps — Figure 6 shows face-on volume density maps. Each map 
shows the H I volume densities in a sheet sliced parallel to the Galactic equator. The interval 
of each map is 500 pc; z = +1000 pc, +500 pc, pc, —500 pc, and —1000 pc. 

3. edge-on volume density maps — Figure 7 shows edge-on volume density maps. Each map 
shows the H I volume densities in a sheet through the Galactic center perpendicular to the 
Galactic plane. Individual panels show sections through 9 = 50° - 230°, 9 = 80° - 260°, 9 = 110° 
- 290°, and 9 = 140° - 320°. 

4-1. Hi Spiral Arms 

Some H I spiral arms can be traced in the resultant face-on maps, figure 5 and figure 
6. In order to investigate the properties of the spiral arms, the H I column density is shown in 
polar coordinates in figure 8. The horizontal and vertical axes represent the azimuthal angle 9 
and the distance R from the Galactic center in logarithmic scale, respectively. Each arm on the 



Z(R,0) =2.27n Hh z 1/2 . 



(13) 
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polar diagram is found to be aligned in a straight line, indicating that the arms are logarithmic 
spiral arms. Schematic tracers of H I arms are presented in figure 9. 

The largest arm is the so-called Outer arm which is traced from (R, 9) ~ (7.0 kpc, — 
125°) to (R,6) ~ (12.1 kpc, 130°). The pitch angle is estimated to be ~ 7°. Another prominent 
arm is traced from (R, 9) ~ (5.6 kpc, 100°) to (R, 9) ~ (11.4 kpc, 310°), which is the so-called 
Sagittarius-Carina arm. Its pitch angle is also estimated to be ~ 11°. Another arm can be 
traced from (R, 9) ~ (5.8 kpc, 50°) to (R, 9) ~ (13.6 kpc, 200°), which corresponds to the so- 
called Perseus arm. The pitch angle is ~ 18°. This arm may apparently connect to the Outer 
arm (figure 5), but has a clearly different pitch angle. Hence, they may be distinct arms. 

The Sun is located at a relatively H i-rich region between the Sagittarius-Carina arm 
and Perseus arm, which is called the local arm. 

The Galactic spiral structure has been studied by many researchers using various meth- 
ods. Oort, Kerr, and Westerhout (1958) found three H I spiral arms: the Sagittarius, Orion, 
and Perseus arms. Although the Orion arm may correspond to the local arm, the others are in 
good agreement with our traced arms. Weaver (1970) found the Outer arm, which also agrees 
with our result. They estimated the pitch angle of the Sagittarius-Carina arm to be 12.5°, 
which is consistent with our result. Kulkarni, Blitz, and Heiles (1982) reported three spiral 
arms, and estimated the pitch angles to be 22° - 25°, which are pretty large compared to our 
values. Georgelin and Georgelin (1976) presented four spiral arms using the Hll regions. They 
traced the Sagittarius-Carina arm with a pitch angle 10° - 15°, and the Perseus arm with s 
pitch angle of 11°, which are consistent with our results. However, their Scutum-Crux and 
Norma arms could not be confirmed in our H I map. Drimmel and Spergel (2001) fitted the 
far and near infrared data form COBE/DIRBE by a 2- and 4- armed logarithmic spiral model. 
Our Sagittarius and Perseus arms are consistent with their spiral structure. Also, our result 
that the Galactic arms are logarithmic spiral arms confirms their logarithmic spiral model. 

4-2. Size and Mass of the Hi Disk 

The radial distribution of the surface density of the H I gas is shown in figure 10 by 
averaging it azimuthally. The surface density is low near the Galactic center, and is lower than 
1.9 M pc -2 within a radius of 6 kpc. The H I surface density attains the maximum values 
at radii 7-12 kpc, where it becomes larger than 4 M Q pc -2 . Then, the H I density decreases 
toward the outer region. It amounts to ~ 1 A'/qPc" 2 at radius 17 kpc, and therefore the size 
of the H I disk is 17 kpc according to the definition of Broeils and van Woerden (1994). This 
radius of the H I disk is about 1.3-times as large as the stellar disk, since the stellar disk is 
estimated to be 13 kpc assuming R = 8 kpc (Robin et al. 1992). This ratio is typical for general 
spiral galaxies. 

The H I mass within a radius of 17 kpc was calculated to be 2.5 x 10 9 M , which 
corresponds to 1.5% of the total mass predicted by the rotation curve. 
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4-3. Warping and Tilting Hi Disk 

Figure 7 shows that the H I disk is strongly warping at the outer region, which is 
consistent with former studies (Burton, te Lintel Hekkert 1986; Diplas, Savage 1991). The 
amplitude of warping is large in the directions 9 = 80° and 9 = 260°. In the direction 9 = 80°, 
warping starts near R = 12 kpc, and the H I disk is bending up to z = +1.5 kpc height from 
the Galactic plane at R = 16 kpc. Also, in the direction of 9 = 260°, warping also starts at 
R = 12 kpc. Then, the H I disk is bending down to z = — 1 kpc height from the Galactic plane 
at R = 16 kpc, and is flopping back to z — —0.5 kpc height at R ~ 18 kpc. Thus, the Galactic 
warping is found to be asymmetric. Since the size of the stellar disk is thought to be about 13 
kpc, the warping starts at the edge of the stellar disk, which follows one of the rules of warping 
(Binney 1992). Hence, the warping is asymmetric. 

This work displays the tilting of the inner disk (Burton, Liszt 1978; Sanders et al. 1984). 
In figure 11, the centroid of the H I disk zq at R = 5 kpc is plotted versus azimuthal angle 9 for 
an instance, which shows that the H I disk is tilted against the 6 = 0° plane. The amplitude 
of the tilting is ~ 100 pc. The H I disk is bending up in the direction 9 ~ 80°, and is bending 
down in the direction 9 ~ 260°, which is the same trend as the warping of the outer disk. 

4-4- Lopsided H I Disk 

Figure 5 shows that the Galactic H I disk is not symmetric. The H I disk is swelled in 
the direction 9 = 310° in the third quadrant. The radius of the outer H I disk (defined by the 
contours at 8.0 x 10 19 H cm -2 in the second quadrant) is ~ 11 kpc, whereas it is ~ 22 kpc in 
the fourth quadrant. Hence, the Galaxy is significantly lopsided in H I. 

4-5. Scale-Height of the Hi Disk 

We define the scale- height of the H I disk as the value of Z\ii (see subsection 3.3). The 
azimuthally averaged scale-height is plotted in figure 4 versus the radius. Figure 4 shows that 
the scale-height increases with radius, from ~ 100 pc at R = 1 kpc to ~ 700 pc at R = 17 kpc. 
According to the Spitzer (1942)'s model, the scale-height is calculated as 



Assuming that velocity dispersion is constant at any radius, the increase in the scale-height 
implies a decrease in the surface mass density at the Galactic plane. 

In subsection 3.1, the volume density of H I gas (n H io) and scale-height (21/2) are obtained 
for each point of the outer Galaxy. In figure 12, we plot the scale-height, Z1/2, versus H I volume 
density, nmo, in an annulus of a constant radius of R = 12.0 ±0.5 kpc. Figure 12 shows that 
there is a negative correlation between the H I volume density, nni , and the scale-height, zu 2 , 
which were determined independently of each other in our analysis. The H I volume density is 
approximately expressed by nmo = 0.63exp(— zi/ 2 /340 pc) by a least-squares fitting. 




(14) 
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4-6. Possible Systematic Problems 

In this work, resultant maps were obtained based on some assumptions, which might 
potentially cause systematic errors in our results. Possible systematic problems are described 
below. 

1. Pure circular rotation — Though the distances of H I clouds were calculated based on an 
assumption that the Galaxy rotates in pure circular rotation, Burton (1972) mentioned that the 
streaming motion predicted by the density wave theory is not negligible. However, because we 
are interested in the global distribution of the H I gas, the effects due to the streaming motion 
were not taken into account. Since the streaming amplitudes were estimated to be between 3 
and 8 km s -1 , the typical errors in the calculated distances of H I clouds were on the order of 
1 kpc. Particularly, in the central region, the non-circular motion may be large, because the 
Galaxy is believed to have a bar whose size is estimated to be about 3 kpc. Hence, the H I 
distribution in the central 3 kpc may be incredible. 

2. Dependence on rotation curves — Our resultant H I map is dependent on the rotation curves 
and there occur two possible systematic errors. First is the H I disk size. If the rotation velocity 
becomes larger values than the one we used, the H I disk will be stretched outwards. Second is 
the H I contrast. If a rotation curve has multiple peaks such as those found in rotation curves 
of Honma and Sofue (1997), H I is concentrated in rings at corresponding radii, because the 
H I density depends on \dV T /dr\ in equation (5). We used a smooth rotation curve in order to 
avoid such artificial rings. 

3. Cylindrical rotation — We assumed that the gas away from the Galactic plane rotates at the 
same speed as the gas in the underlying disk. This assumption was often used in former studies 
(Lockman 1984; Diplas, Savage 1991). However, this assumption may not be adopted for the 
high-altitude H I gas, because it is known that the rotation velocity at a large distance from 
the Galactic plane becomes slower than that of the gas at underlying disk (Swaters et al. 1997) 
based on observations of an external edge-on galaxy, such as NGC 891 in H I line. 

4. Fitting method — We determined the H I distribution in the inner Galaxy by a fitting 
method, as mentioned in subsection 3.3. However, this fitting method cannot always work 
precisely. Largely deviated H I clouds from the surrounding gas, often seen in resultant maps, 
are thus doubtful, while a global feature is plausible. 

5. Dependence of adopted functions for z- distribution — While the function sech 2 (2;) was 
adopted in our analysis, the real gas distribution can not be perfectly represented by it, because 
the gas is not isothermal (Lockman 1984). We present H I column density maps derived by using 
other functions, such as a Gaussian function and exponential function in figure 12. Though the 
resultant maps depend on the adopted functions, the global structure is not strongly affected 
by the difference of the chosen functions, but the resultant densities varies by a factor of ~ 2 
inside the solar circle. 
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5. Summary 



We obtained a three-dimensional H I distribution in the Milky Way Galaxy using the 
Leiden-Dwingeloo H I survey (Hartmann, Burton 1997), Parkes survey (Kerr et al. 1986), and 
the NRAO survey (Burton, Liszt 1983) by combining with Clemens (1985) and Dehnen and 
Binney (1998) 's rotation curves. The distance and density of H I clouds were determined by 
assuming that the Galaxy rotates circularly. We solved the near-far problem in the inner 
Galaxy using a fitting method based on a model of the vertical distribution (Spitzer 1942). 

Our resultant map shows: 

1. The Galaxy has three prominent spiral arms: the Outer arm, the Sagittarius-Carina arm, 
and the Perseus arm. All of the arms are logarithmic spiral arms. 

2. The H I disk is strongly warping in the outer disk and the warping is asymmetric. Also, the 
inner disk is found to be tilting. 

3. The size of the H I disk is about 17 kpc, and the H I mass within this radius is about 2.5 x 10 9 
M . 

4. The H I disk is swelling in the southern hemisphere, which indicates that the Galaxy is 
significantly lopsided. 

5. The scale-height increases with the radius globally. The scale-height is correlated with the 
H I volume density at the centroid of the H I layer. 
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Figure Captions 



Fig.l. Rotation curve used in this paper. The rotational velocity in range of R < 8 kpc is cited 
from Clemens (1985), and that in range of R > 8 kpc is cited from Dehnen and Binney (1998). 

Fig. 2. Radial velocity versus the heliocentric distance in the direction of / = 30° and b = 0°. 

Fig. 3. Radial velocity field of the Galactic plane calculated with the rotation curve. The solid 
lines denote positive velocities and the dashed lines denote negative velocities. 

Fig. 4. Scale-height versus the Galactocentric distance. 

Fig. 5. H I column density map. The contours are drawn at levels of 8.0e+19, 1.6e+20, 3.2e+20, 
6.4e+20, 1.3e+21 atom cm" 2 . 

Fig. 6. Horizontal cross-section maps. Each map shows the H I volume densities in a sheet 
sliced parallel to the Galactic equator. The interval of each map is 500 pc; z = +1000 pc, +500 
pc, pc, —500 pc, and —1000 pc. The contours are drawn at levels of 0.01, 0.02, 0.04, 0.08, 
0.16, 0.32, 0.64 atom cm" 3 . 

Fig. 7. Vertical cross-section maps. Each map shows the H I volume densities in a sheet through 
the Galactic center perpendicular to the Galactic plane. Individual panel shows sections through 
9 = 50° - 230°, 9 = 80° - 260°, 9 = 110° - 290°, and 9 = 140° - 320°. The contours are drawn 
at levels of 0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 0.64 atom cm" 3 . 

Fig. 8. H I column density map in polar coordinates. The vertical axis denotes the 
Galactocentric radius on the logarithmic scale, and the horizontal axis denotes the azimuthal 
angle around the Galactic center. The contours are drawn at levels of 1.0e+20, 2.0e+20, 
4.0e+20,8.0e+20, 1.6e+21 atom cm" 2 . 

Fig. 9. (a)Schematic tracers of H I arms superposed on a resultant H I column density map. (b) 
Schematic tracers on a polar-coordinates diagram. 

Fig. 10. Surface density of the H I gas versus the Galactocentric distance. 

Fig. 11. Plots of centroid zq of the H I disk versus the azimuthal angle 9 around the Galactic 
center at R = 5 kpc. 

Fig. 12. Changes in the scale-height at an annulus of constant radius R = 12.0 ±0.5 kpc. The 
horizontal axis denotes the H I density at the centroid of the H I layer (nHi ) and the vertical 
axis denotes the scale-height (21/2)- 
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Fig. 13. (a) H I column density map derived by using a Gaussian function to determine the 
H I distribution in the inner Galaxy, (b) H I column density map derived with an exponential 
function. The contours are drawn at levels of 8.0e+19, 1.6e+20, 3.2e+20, 6.4e+20, 1.3e+21 
atom cm -2 . 
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